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A one-dimensional model with water-like anomalies and two phase transitions 
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We investigate a one-dimensional model that shows several properties of water. The model 



Q-i combines the long-range attraction of the van der Waals model with the nearest-neighbor 

■ interaction potential by Ben-Nairn, which is a step potential that includes a hard core and 

. a potential well. Starting from the analytical expression for the partition function, we 



determine numerically the Gibbs energy and other thermodynamic quantities. The model 
shows two phase transitions, which can be interpreted as the liquid-gas transition and a 



S transition between a high-density and a low-density liquid. At zero temperature, the low- 

density liquid goes into the crystalline phase. Furthermore, we find several anomalies 
that are considered characteristic for water. We explore a wide range of pressure and 
temperature values and the dependence of the results on the depth and width of the potential 
7^3 ■ well. 
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I. INTRODUCTION tulated 20 years agcA Despite the fact that the 

idea of characterizing liquid water by a high- 

rv, 1 Due to its omnipresence on earth and its im- , , . , . , , 

r density and a low-density liquid is old and has 

go : portance for life, understanding the behavior of ^ ^ this so . called LDL-HDL 

water and water mixtures is of utmost impor- . . ...... . . 

Q\ 1 r transition, which implies that water has a sec- 

2 ' tance. Although the anomalous properties of , . . , . . .„ . , „ 

cN ■ ond critical point, is still controversial. Sev- 

water have been subject to investigation for a , lA , , , 

J ° eral alternative scenarios have been suggested, 

long time, the reasons and underlying princi- r , . , , , ...... . , 

^ ■ J ° r all oi which can be described within a simpl 

pies leading to the outstanding characteristics of 

the model parameters 5 . Nevertheless, a LDL- 
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v ° ° cooperative hydrogen-bond model by varying 

water are not yet completely understood 1 ' 2 . Re- , , , 5 



cent effort has focused on the properties at low TT _ T ... . . _ 

HDL transition is a generic scenario for core- 
temperatures and high pressures, where many ^ , . , 6 _ 8 , . . ^ „ 
r softened potentials , and the experimentally 

different ice phases were found, and where a , , , • , , r 

observed dynamical crossover between fragile 

liquid-liquid phase transition between a high- , .... , . . 

n n r and strong behavior in supercooled water can be 

density phase and a low-density phase was pos- . , c .. „ T ., Q 

J r J r r interpreted as a signature of the Widom hne\ 

a) lotta@fkp.tu-darmstadt.de „. , ^ _ „ , 

Theoretical studies of water range from be- 



ing as precise as possible to being as simple as 
possible. One older example for a description 
of the behaviour of water that allows for the cal- 
culation (and thus prediction) of properties with 
large precision is an equation of state for wa- 
ter with 58 parameters 10 . With such complex 
descriptions, however, no mechanisms can be 
identified, and the underlying principles leading 
to the characteristic behaviour of water remain 
unrevealed. Nowadays, the dominant theoreti- 
cal approaches to water are molecular dynam- 
ics simulations and Monte Carlo simulations, 
which have become possible due to increasing 
computing power, allowing for a microscopic 
modelling of water starting from the full quan- 
tum mechanical description of its molecules, or 
from simpler, classical molecular models with 
few parameters. These models were able to re- 
produce many properties of water—. 

However, critical voices 12 argue that com- 
puter simulations of water do not follow the 
basic guidelines of theoretical investigations. 
Though promising in many respects, the sys- 
tems consist of hundreds to thousands particles 
interacting via semi-empirical interaction po- 
tentials, for which the interpretation and iden- 
tification of mechanisms remains difficult. A 
need has been identified to find even simpler 
models capturing relevant features of water in 
order to understand which ingredients are nec- 
essary and indispensable for the occurrence of 
anomalies or other features of real water. 



An important step in the understanding of 
water was the insight that neither the bent shape 
of the molecule nor the directedness of hydro- 
gen bonds are necessary for many of the ob- 
served features of water. In fact, a soft core 
in a radially symmetric interaction potential is 
sufficient to allow for a discontinuous change 
in the preferred distance between molecules as 
the pressure is varied, leading to various water 
anomalies 13 . The idea of modelling water with 
an effective soft-core potential is experimentally 
supported by potentials derived from the exper- 
imental 0-0 radial distribution function-^. 

A further simplification is implemented by 
the Jagla potential, which has no attractive 
part at all, but a ramp that is a simple repre- 
sentation of a soft core. Even simpler are re- 
pulsive step potentials which were introduced 
in the 1970's^ and are still investigated un- 
til today-^— , although there are results show- 
ing that a single temperature- and pressure- 
independent potential may not be sufficient to 
capture all properties of water—. 

Even one-dimensional water models have 
been used to describe water anomalies. A 
lattice-gas model, where distance 2 between 
nearest neighbors is associated with a stronger 
binding energy than distance 1, shows a zero- 
temperature critical point and thermodynamic 
and dynamic anomalies similar to water—. A 
rather complex lattice model with two differ- 
ent repulsion scales and a mean field attraction 
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leads to a rather rich phase diagram— Con- 
tinuous one-dimensional models with a short- 
range soft-core potential produce various wa- 
ter anomalies 24 , which occur also in the two- 
dimensional version 2 ^. A model with two wells, 
which was first published in 1992^, was ar- 
gued to provide an explanation for the den- 
sity anomaly^ 7 -. The simplest one-dimensional 
model that produces water anomalies was in- 
troduced by Bell— and re-discovered by Ben- 
Nairn—. It contains a step potential with a sin- 
gle well and shows among other properties a 
density anomaly and a minimum in the isother- 
mal compressibility. A drawback of such one- 
dimensional models with short-range interac- 
tions is, however, that they cannot show a phase 
transition at temperatures larger than zero^. 
Several general reviews about simple and sim- 
plified water models have been written 2 *^-, and 
a very good introduction into water models can 
be found in the book by Ben-Nairn 3 ^. 

In addition to these newer models for wa- 
ter, there exists also the phenomenological mean 
field approach of van der Waals, introduced in 
his PhD thesis in 1873, that has since then been 
discussed extensively. The properties of a one- 
dimensional model in the van der Waals limit, 
where the range of the force goes to infinity 
while its strength goes to zero, were discussed 
in the 1960's by several authors^— . Ideas for 
adjusting this model to give a better description 
of water were also discussed by Heidemann and 



PrausnitzaS, where a van der Waals model for 
fluids with associating molecules is introduced. 
It was noted, however, the van der Waals model 
generally fails for water because there are strong 
directed interactions in water that can not be 
modeled with a mean field attraction 37 . 



In this article, we will study a one- 
dimensional model that combines the features 
of the Ben-Nairn model and the van der Waals 
model. Combining a simple water model with 
a long-range attraction was first done by EA. 
Jagla 3 ^. He found that by adding a long-range 
attraction similar to the van der Waals gas to 
his above-mentioned three-dimensional model, 
the water-like properties are maintained and an 
additional phase transition can occur. Adding 
a long-range attraction to a one-dimensional 
model, as will be done in this paper, leads to 
the probably simplest possible one-dimensional 
model that shows two phase transitions and vari- 
ous water anomalies. The mean-field term turns 
the zero-temperature phase transition of the 
Ben-Nairn model, which is due to the potential 
well, into a finite-temperature phase transition 
between two different liquid phases. Further- 
more, the mean-field term introduces the liquid- 
vapor phase transition into the one-dimensional 
model. 
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FIG. 1. The short-range interaction potential for a 
ID model of water suggested by Ben-Nairn^ 9 -. 

II. MODEL 

Our model is a combination of the one- 
dimensional model by Ben-Nairn and the van 
der Waals model. We first describe these two 
ingredients of the model separately before we 
present the combined model. 

A. Ben-Nairn model 

The Ben-Nairn model was first analyzed by 
Bell in 1969^ and discussed later with respect 
to water by Ben-Nairn 29 . This model describes 
a one-dimensional system of N particles inter- 
acting via a short-range potential. The potential 
sketched in figure Q] consists of a hard-core re- 
pulsion (corresponding to the excluded volume 
of a particle) and a minimum that mimics the 
effect of hydrogen bonding. The corresponding 
partition function can be calculated analytically, 



gBN(T, p) 
Iln. 



(1) 



fipA 



j3 Q-Ppo-i + Q-Bpo- 2 (Q-Be - 1)(1 - Q-PPS) 

and to the specific volume 

v B n(T,p) = ^-\ T (2) 
dp 

with /3 = j and the thermal wavelength A = 
h/ ^InrnksT. This model exhibits a density 
anomaly and a minimum of the isothermal com- 
pressibility, in agreement with two character- 
istic features of real water. Additionally, for 
low temperatures a steep transition from a lower 
density to a higher density is observed with in- 
creasing pressure. At T = 0, this transition be- 
comes a real first-order phase transition, with all 
particles having distance <j\ for p(cr 2 - en) > \e\ 
and distance <x 2 otherwise. This phase transition 
can be calculated by determining the phase that 
minimizes G = E - TS + pV at T = and for 
the considered value of the pressure p. 

B. Van der Waals model 

The van der Waals gas is described by the 
thermal equation of state 

(p + ^)(v-b) = k B T, (3) 

where a and b are parameters of the considered 
gas and v is the specific volume v = V/N. For 
our model, we have b = cr\. The pure van der 
Waals model leads to the well-known phase dia- 
gram of the van der Waals gas with a first order 



phase transition line between a fluid and a gas 
terminating in a critical point, which is situated 
at T c = $a/27k B b, p c = a/llb 2 and v c = 3b. If 
a ± and b + 0, the thermal equation of state is 
equivalent to 

p + ij(3v-l) = 8f (4) 

with the reduced variables v = v/v c , f = 
T/T c , p = p/Pc chosen such that the critical 
point isatv = j5 = 7' = l generically^. Al- 
though this model has a liquid-gas phase transi- 
tion, no other features of water are matched. 

C. The combined Ben-Naim-van der 
Waals model 

We now combine the two ingredients and in- 
vestigate the resulting model, which we call the 
Ben-Naim-van der Waals model. Below, we 
will use the index BNJ for this model, because 
Jagla first introduced a van der Waals term into 
a water model. The short-range interaction po- 
tential has the same form as the potential used 
by Ben-Nairn (see figure [T]). Additionally, we 
introduce a long-range interaction analogously 
to the van der Waals model via the substitution 
p — » (p + a/v 2 ) in the above expression for the 
specific volume of the Ben-Nairn mode l 15i4 ° i41 . 
The reduction of pressure by a/v 2 follows di- 
rectly from a decrease of the internal energy by 
a/v due to the long-range attraction. We there- 



fore obtain 

a 

v(T, p) = v B nj(T, p) = v BN (T, P + -?) (5) 

v 1 

a la 

g(T,p) = g B Nj(T,p) = g BN (T,p + — ) . 

V 1 V 

(6) 

A more detailed derivation of expression © 
is given in Appendix B in the publication by 
Truskett and Dill 4 ^. The implicit equation for 
the specific volume can be solved numerically. 
We use again the reduced variables (p, f, v) and 
introduce further the dimensionless variables 

Cfj = crJVc , CT 2 = <T 2 /V C , 5 = 6/v c , £ = 

e/k B T c , g = g/k B T c , which are chosen such 
that lengths are measured in units of v c and 
that for e = 0, the model reduces to the pure 
van der Waals system with a critical point at 
v = T = p = 1 . Note that since T c oc a, we have 
6 oc e/a and the reduced variable e measures 
thus the ratio between the depth of the short- 
range potential well e and the global attraction 
strength a. 

This model has now three remaining free pa- 
rameters e, cf 2 and 6. Note that since the hard 
core distance is b = cr l5 we have &\ = 1 /3, and 
this parameter has thus been eliminated. From 
now on, we will omit the tilde, and we will set 
x = x for the variables x = T,p, v, cr\ , cr 2 , 6, e, g, 
which is equivalent to using dimensionless vari- 
ables as defined above. This choice is not pos- 
sible for the pure Ben-Nairn model, where other 
(arbitrary) units are used. 
The expression for the specific volume v(T,p) 



can be evaluated numerically in order to obtain 
a density profile. Whereever there is more than 
one solution for v, one must determine the solu- 
tion that minimizes the Gibbs energy in order to 
obtain a phase diagram. We also evaluated the 
isothermal compressibility and the isobaric heat 
capacity 

1 dv 

kt = — T-lr (7) 
v op 

ds ^§ 

The first can be extracted from v(T, p), and the 
second follows from equation ©. Additionally, 
the probability density of the distance between 
nearest neighbours can be calculated. It is iden- 
tical to the weight occuring in the expression for 
the partition function and is given by 

U(r) 3(p+3/v 2 ) 

q(r) = ae T e 8r (9) 
where or is a normalization constant ensuring 

J~oo 
q(r)dr = 1. 

IH. RESULTS 

If not noted otherwise, the parameters of the 
model are cr 2 = 2/3 and 8 = 1/10. These values 
are chosen such that their ratio matches that of 
the parameters used by Ben-Nairn 29 (while we 
have to set cr\ = 1 /3). The parameter e is varied 
from down to -10. 



A. Density profiles and phase diagrams 

Figure [2] shows the density profile, the co- 
existence regions of two phases, and the phase 
diagram in the p-T plane for different values of 
6. The bottom row shows the pure Ben-Nairn 
model with cr, = 1/3, cr 2 = 2/3, 6 = 1/10 
and £ = -5. For small pressures and temper- 
atures, the specific volume is v « 0.67, imply- 
ing that most nearest-neighbor pairs are hydro- 
gen bonded. At high pressure, the specific vol- 
ume decreases to v ~ 0.34, where the nearest- 
neighbor distances are close to the hard core di- 
ameter. At T = 0, this transition is a real phase 
transition and occurs at p = -e/(cr 2 - <T\) = 15. 
The top row of Figure 1 shows the pure van der 
Waals model (where e = 0). There is a first or- 
der phase transition (VdW-transition) between 
a phase with a small volume (corresponding to 
a high density fluid phase) and a phase with 
a higher volume corresponding to a low den- 
sity gaseous phase. As temperature decreases 
or pressure increases, the density of the liq- 
uid phase increases continuously, however with- 
out displaying a region of particularly steep in- 
crease, as is the case in the Ben-Nairn model. 

The rows between the first and last show the 
changes in the system as |e| ist increased. The 
first change that can be seen is the appearance 
of a metastable phase at very low temperatures. 
This phase has a higher density than the stable 
phase, which is influenced by the presence of 
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FIG. 2. Density profiles (left), coexistence regions (center), and phase diagrams (right) for e = 
0,-1, -3, -5, -10, -oo (from top to bottom). The color code in the density profiles (left column) shows 
the specific volume v for each state point, while in the coexistence region (central column) indicates the 
absolute value of the difference |Ag| in the Gibbs energy between the coexisting phases. In the bottom line, 
which represents the pure Ben-Naim model, unsealed parameters p and T are shown, and the value e = -5 
was used. Since the Ben-Naim model shows no phase transition, the location of steepest change in the 
specific volume in the "phase diagram" is indicated by a dotted line. 
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the potential well. With increasing |e|, the re- e = -3, where both phase transitions are pro- 
gion where the metastable phase exists becomes nounced. Figure |4] shows the distribution q(r) 
larger (see Figure [3] for a zoom into the |e| in- (see Equation ©) for different temperatures 
terval between -1 and -2), and at e ~ -1.5 and pressures corresponding to states below and 
the metastable phase becomes stable for the first above the LDL-HDL-phase transition. The dis- 
time. This first occurrence of a phase transition tribution broadens with increasing temperature, 
between the two liquids can best be understood and the probabilities for two particles to have 
by considering the case T = 0: At T = and the H-binding distance or the smaller distance, 
going back to the original units, the Gibbs en- respectively, become more similar. Figure \5\ 
ergy per particle in the high density phase is shows the proportion of nearest neighbor pairs 
gh = per i -a/cTi, and it is gi = e + pcr 2 -a/cr 2 in in each of the three distance intervals as a func- 
the low density phase. This corresponds to g h = tion of T and p. One can clearly see that for 
p(Ti/S -27/80"! and gi = e + pcr 2 /& -27/8cr 2 in high temperatures and low pressures, most par- 
reduced variables. Since pcr 2 > pcr u the Gibbs tides have a distance larger than cr 2 + 6, while 
energy of the high density phase can become for low temperature and pressure most particles 
smaller than that of the low-density phase only sit in the potential well, and for low T and high 
when p most particles have the minimum distance <j\ . 

27/1 i \ We end this section by discussing in more de- 

-e>— = 1.6875 (10) 

8 \ct\ cr 2 ) tail the LDL-HDL transition of our model. The 

This is why the second phase transition occurs slo P e of the LDL-HDL phase transition should 
only for sufficiently large absolute values of e have a negative sign for water, since the low- 
and extends to T = only if 6 < -1.6875. density phase is associated with a higher degree 
When e increases further, the second transition of order ' while the high-density phase is less or- 
is shifted towards higher pressures, and for e = dered - This is different from most other mate- 
-10 the transition is already out of the shown rials and is rather difficult t0 obtain in a sim P H - 
parameter range. The first (VdW-) transition fied mode1 ' because in two or more dimensions 
is shifted towards higher temperatures with in- a lar § er avera § e P article distance corresponds to 
creasing e more available states. This is not true, however, 
In order to better visualize the properties of for our one-dimensional model. The slope of 
the system on both sides of the phase transitions, the P hase transition can easily be changed in this 
we show in Figures @ and the equilibrium model b y adjusting the width of the well, 6, as 
distance distribution of particles for the value is shown in Fi § ure ® The slo P e of the P h ase 
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FIG. 3. Coexistence regions for the combined Ben-Naim-van der Waals potential for —e - 1 , 1 .25, 1 .5, 1 .75 
(from (a) to (d)). 

B. Anomalies 



The combination of a hard core and a 
potential well at a larger distance leads to 
various anomalies similar to those of water, 
as was shown by Ben-Nairn. The number of 
anomalies of water is large, see for example 
http : //www . btinternet . com/~martin. chaplin/anmlie 
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and not all of them can occur in a simple model. 
FIG. 4. Probability q(r) for a distance r between two In tne following, we investigate three different 
neighbouring particles for e = -3.0 and different val- anomalies in our model, which are also present 
ues of T and p. in the original Ben-Nairn Model. We choose 

again e = -3.0, since for this value both phase 
transitions are well visible, 
transition is negative for sufficiently small 8. The density anomaly is shown in figure [71 



(a) 




(O 

FIG. 5. Probability q{r) for e = -3.0 and for a 
distance r e [cri,o" 2 ] (a), r e [o~2,cr 2 + 5] (b), 
r e [cr 2 + 8, oo] (c). 

The specific volume has a region of negative 
slope and a local minimum, where density is 
maximum. The region of negative slope, where 
the thermal expansion coefficient 

a P = -^\ P (U) 
v oT 

is negative, becomes broader with increasing p, 
due to the closeness of the second phase transi- 
tion. The right boundary of the red area cor- 
responds to a local density maximum, while 
the left boundary corresponds to a local den- 
sity minimum (as function of T). With increas- 
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FIG. 6. Density profiles for e = -3.0 and (a) 6 = 
0.01 and (b) 6 - 0.002. 

ing temperature, more and more particles move 
out of the potential well and preferentially to its 
left side (see also figure [5]). Beyond the red re- 
gion, more and more particles move to the right 
side. Above a pressure of p ~ 10, the den- 
sity anomaly vanishes. This behaviour is also 
known from real water, which does not show 
anomalies at high pressures 43 . 

The isothermal compressibility and the iso- 
baric heat capacity are shown in Figure [H For 
a simple hard-rod system, one would expect 
that kt(T) is a monotonically increasing func- 
tion, while kt(p) is a monotonically decreas- 
ing function. As to c p , this would be a con- 
stant for all pressures and temperatures in a hard 




figure |8(a)| and |8(b)[ ). The maximum for high 
pressure (p = 20) at T « 0.8, which is visible 
in both quantities, is due to the closeness of the 
Ben-Nairn transition. 

The maxima in K T (p) and c p (p) (see figure 
8(c)| and |8(d)[ ) are also due to the Ben-Nairn 
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FIG. 7. Density anomaly for e = -3.0. Figure (a) 
shows v(T) exemplarily for p = 0.64 and (b) shows 
the region where the coefficient of thermal expansion 

a p = v§f\p is negative- 
rod system^. In our model, due to the pres- 
ence of the potential well, both quantities show 
as a function of T a minimum at low pressure, 
which vanishes at high pressure, as for the den- 
sity anomaly. The singularities in the curves 
Kj(T) and c p (T) for low pressure correspond to 
the van der Waals transition, they are shifted to 
higher temperature for increasing pressure and 
turned into maxima above the critical point (see 



transition, as can be seen by comparison with 
the phase diagram. For the high temperature 
T = 1.5, the van der Waals transition is crossed, 
leading to singularities at p « 2. 

The features corresponding to the Ben-Nairn 
transition were observed similarly in the pure 
Ben-Nairn model^. For instance, the minima of 
kt(T) and c p (T) for small temperatures are also 
present without the van der Waals term. The 
maximum of K T (p) at the BN-transition that was 
observed in the pure model is still present above 
the critical point at the so-called Widom line, 
but shows a discontinuity in the region of the 
phase transition as is expected for a real phase 
transition. The modulation of c p (p) that was 
observed for the pure model^ around the BN- 
transition is not visible in our data, but this may 
be due to the lower resolution we employ for 
the pressure. The overall behaviour of c p (p) 
shows also the signatures of both phase tran- 
sitions. The signatures of the van der Waals 
transition (discontinuities and maxima in the re- 
sponse functions) were, of course, not seen in 
the pure model. 

To illustrate the behaviour of the response 
functions k t and c p more globally, we show in 
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FIG. 8. Left column: Isothermal compressibility kj as a function of T (a), as a function of p (c,) and regions 

where kt as a function of T has a negative slope (e). Right column: Isobaric heat capacity c p as a function 

of T (b) and as a function of p (d), and regions where c p as a function of T has a negative slope (f). All 

plots are for a potential with e = -3.0. 



12 



figures |8(e)| and |8(f)| the regions where these sity has a similar shape as the curve we find for 



variables have a negative slope as functions temperatures above the BN critical point, 
of the temperature T. Just as for the den- The anomalies of our simple one- 

sity anomaly, negative slopes dx T /dT < and dimensional model resemble thus in several 

dc p /dT < are considered as 'anomalous' be- respects those of real water and water models 

haviour. One can see that both functions have although not all properties of water anomalies 

three regions of negative slope: The region at are matched. We want to emphasize two 

small T and small p, where the potential well features that seem to be particularly interesting 

traps particles that would otherwise be to its left; with respect to real water. First, the region of 

the region close to the HDL-LDL phase transi- anomaly of k t (T) and c p (T) is restricted to low 

tion and around the Widom line; and third the pressures and the system behaves normally at 

region above the van der Waals liquid-gas tran- high pressures, which is similar in real water- 4 ^. 

sition. The negative slopes next to the phase Second, the anomalies are not independent 

transitions are expected since both k t (T) and of the phase transitions. We observed that 

c p (T) have a maximum at the phase transitions, the coefficient of thermal expansion a p (T) is 

A classification of parameter regions of strictl y positive as long as e is close enough 

anomalous properties of real water was made to zero such that there is no BN-transition. 

by Errington and Debenedetti 44 , who also in- However, when a phase transition occurs (as for 

vestigated the interdependence of structural, dy- e = ~ 3 -°> for example), the density anomaly 

namic and thermodynamic anomalies. The re- is also seen in a temperature region above the 

gion of thermodynamic anomalies, where the critical P oint of this transition (compare figure 

density increases upon increasing temperature, El and [7(b)]). This supports the hypothesis that 

spreads over a smaller density interval when anomalies of water may be related to phase 

temperature increases. Correspondingly, we transitions at a lower temperature that is not 

find in our model that with increasing temper- observable due to spontaneous crystallization, 

ature, the region of pressure where an anomaly but S ives rise to signatures in form of anomalies 

occurs gets narrower. The region of anomalous at ambient conditions, 
density behaviour of three-dimensional parti- 

cles interacting via isotropic core-softend po- Iy DISCUSSION AND CONCLUSION 

tentials was investigated by de Oliveira and 

coworkers 45 . The line in the p-T diagram corre- We have introduced a simple one- 

sponding to the temperature of maximum den- dimensional model that combines two existing 
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models and shows several properties of real wa- 
ter. The model has the following two features: 
1.) A short-range potential that introduces 
two length scales and correlates low binding 
energy with a low density and thus with an open 
structure^ and 2.) a long-range interaction 
that lowers the energy when the total volume 
is smaller. These are two properties that real 
water possesses: The short range potential 
leads to the formation of hydrogen bonds that 
correlate low binding energy with an open 
structure. Additionally, there is a long-range 
interaction from electrostatic interactions and 
van der Waals forces that encourages high 
density. In contrast to the pure Ben-Nairn 
model, which has only short-range interactions 
and therefore no phase transition, our model has 
two first-order phase transitions, a liquid-gas 
transition and a HDL-LDL transition. An 
important parameter of the model is the ratio 
of the two energy scales that are given by the 
depth of the potential well and by the strength 
of the mean- field attraction. Only when this 
ratio 6 is high enough does the second phase 
transition occur. With increasing e, the two 
phase transitions move to higher T and higher 
p, respectively. The relative importance and 
balance between the long-range van der Waals 
force and a short range hydrogen bonding is 
presently discussed for ice4£. Also in a one- 
dimensional lattice model the balance between 
the mean field attraction and the short-range 

14 



interaction potential had to be adjusted in order 
to obtain a phase diagram similar to water—. 

Of course, a one-dimensional model cannot 
be a complete and realistic description of real 
water, and it cannot make quantitative predic- 
tions. However, a simple model can help to 
gain insight into the principles underlying the 
special properties of water, and it may lead 
to more understanding than a more complex 
model. The partition function and implicit ex- 
pressions for various thermodynamic quantities 
could be given analytically, and were evaluated 
numerically. The results obtained using simple 
models, which can be written down analytically, 
can be interpreted much easier than those ob- 
tained with more complex models, which would 
require computer simulations 31 . 

A generalization of our model to higher di- 
mensions could be done in several ways, all of 
which have their drawbacks. In particular, the 
geometry must be considered. If the potential 
would be generalized to three dimensions by 
making it spherically symmetric, the number of 
states in the potential well would be much larger 
than at shorter distances, leading to entropic ef- 
fects. For example, the slope of the HDL-LDL 
transition would almost certainly be positive be- 
cause of the Clausius-Clapeyron relation. In in- 
vestigations of the two-scale Jagla-Potential in 
three dimensions, for example, a LDL-HDL- 
transition was found, but also with a positive 
slope in the p-T-plane^-. 



In order to compensate for these entropic 
effects, angle-dependent potentials would have 
to be introduced as for example in a very re- 
cent publication by the Stanley group 48 , where a 
Widom line with a negative slope was found for 
strong tetrahedral interactions. This approach, 
however, is already very close to the detailed 
modelling of single molecules, where a tetrahe- 
dral geometry emerges from the presence of two 
hydrogen atoms arranged in a nearly tetrahedral 
angle with respect to oxygen, and thus provides 
not much more insight compared to molecular 
models. Also, it would be necessary to know 
whether the four sites where H-bonds can be 
formed are independent. This seems to be the 
case according to Predota et al.— . 

The relation between the dimension of the 
model, anomalies, phase transitions and the in- 
teraction potential is discussed by Buldyrev et 
al.— . In this article, a double-step-potential is 
investigated in one, two and three dimensions 
and also the difference between two- and three- 
dimensional models is discussed. The authors 
find in their model systems that liquid anoma- 
lies and a liquid-liquid phase transition may oc- 
cur independently and that a density anomaly 
in a low density phase is not seen when the 
LDL-HDL phase transition line has a positive 
slope. We can not confirm these effects as gen- 
eral trends, as we observe the density maximum 
only when a LDL-HDL transition is present and 
it occurs in the LDL phase. However, we agree 
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about the shape of the temperature of maximum 
density T p (p) and about the fact that this func- 
tion has a maximum. 

Models in more than one dimension do 
not require a long-range attraction in order to 
show phase transitions. Therefore, most higher- 
dimensional models have only short-range po- 
tentials, but all of them require two different 
length scales. Indeed, it has been discussed by 
various authors that a two-length- scale potential 
is a necessary ingredient in order to obtain ther- 
modynamic anomalies similar to water— and 
there is evidence that the hierarchy of anoma- 
lies is determined by the relation between the 
two length scales^. The exact form of the two- 
length- scale potential seems not to be impor- 
tant for the occurrence of anomalies, since they 
occur also in a model with a repulsive step^ 
instead of an attractive well as in the present 
manuscript. In any case, the interaction poten- 
tial has to be such that a higher density is corre- 
lated with a higher energy and a lower density 
with a lower energy, and our work confirms that 
this correlation is necessary for the presence of 
anomalies. However, water is a complex liquid 
and there is some evidence that simple princi- 
ples for the interaction potential may not be the 
whole story. Errington et al.— and later Yan et 
al.— suggested that the occurrence of anoma- 
lies can be related to and predicted by the ex- 
cess entropy, implying that this quantity may be 
more relevant for the occurrence of anomalies 



than the shape of the interaction potential. It 
has to be noted, however, that some assumptions 
of Errington's work do not hold in anomalous 
regions 5 -. Another recent article on this topic 
states that a two-length- scale potential may not 
fully account for anomalies, but that energetic 
and entropic effects may be relevant as well 5 ^. 

As a final note, we would like to point out 
that there exist other materials besides water that 
exhibit anomalies. Anomalies in silica are for 
example investigated by Shell et al.— , and in 
the work by Hoye and Lomba 23 the compari- 
son of water with other tetrahedral substances 
such as Si or Ge was made. In a publication 
by Angell et al. 58 , further tetrahedral liquids are 
mentioned. The simple model discussed in this 
paper, as well as other simple models, do thus 
not only help to understand water, but also other 
materials that have similar anomalies. 
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